Supersymmetry, shape invariance and the Legendre equations 
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In three space dimensions, when a physical system possesses spherical symmetry, the dynamical 
equations automatically lead to the Legendre and the associated Legendre equations, with the 
respective orthogonal polynomials as their standard solutions. This is a very general and important 
result and appears in many problems in physics (for example, the multipole expansion etc). We study 
these equations from an operator point of view, much like the harmonic oscillator, and show that 
there is an underlying shape invariance symmetry in these systems responsible for their solubility. 
We bring out various interesting features resulting from this analysis from the shape invariance point 
of view. 
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Supersymmetry has been shown to be an important 
tool for understanding physical systems of current inter- 
est. It is mandatory for superstrings, and it finds ap- 
plications in diverse areas of physics. A very nice ex- 
ample of this concerns the study of topological defects: 
Bogomol'nyi [l| and Prasad and Sommerfield Q have 
shown that under special circumstances, kinks, vortices 
and monopoles appear as stable, minimal energy states 
that solve first-order differential equations, forming BPS 
configurations in a supersymmetric environment. The 
BPS states are of current interest, since they offer differ- 
ent lines of investigations in physics. For instance, in the 
case of kinks in two-dimensional space-time, the study of 
classical stability of BPS states nicely leads us to super- 
symmetric quantum mechanics (SUSYQM) There 
is also the inverse route where we can re-construct BPS 
states starting with a given SUSYQM system, see, e.g., 
Refs.fl-il. 

For conventional Hermitian Hamiltonians (of the type 
H =p'^ + V{x) with 2m = 1 = in one dimension), the 
SUSYQM problem is described by the first-order opera- 
tors 
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Here the real function W{x) is known as the superpo- 
tential, and with the above operators, we can define two 
Hermitian Hamiltonians 



dx^ 
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+ V-{x), 



= AA^ = -^ + V+{x), 



V^{x) = W^x) T 



dW{x) 
dx 



(2) 



These two Hamiltonians are supersymmetric partner 
Hamiltonians and share all the eigenvalues except for the 
ground state of which is assumed to have a vanishing 
eigenvalue (almost isospectral). The lowering and rais- 
ing operators. A, A'' , also take eigenstates of to those 



of ff+ and vice versa. Such a construction can also be 
carried out for non- Hermitian Hamiltonians 

Supersymmetry relates only a pair of Hamiltonians and 
says that they are almost isospectral without actually de- 
terminingthe spectrum. On the other hand, shape in- 
variance 1^, leads to a sequence of Hamiltonians which 
are pairwise supersymmetric. If the potential depends on 
a real parameter (e.g. a coupling constant) in such a way 
that, 



V+{x,ao) = V-{x,ai) + R{ai), 



(3) 



where oi = /(ao) is a given function of the parameter oq 
and R{ai), representing the spacing between the ground 
states of the two partner Hamiltonians, is a known (con- 
stant) function of the parameter, the potentials (Hamil- 
tonians) are said to be shape invariant. In this case, we 
can construct a sequence of Hamiltonians [IQ] 

i7(o)=i/_(ao), 

H'-'^ = H+{ao) = H^iai) + Riai), 

^ i7_(a2) + R{ai) + R{a2), , 

k 

ffW =i7_(afc) + ^i?(aj),-- - , 



(4) 



all of which will be pairwise supersymmetric as shown in 
Fig. [T] Shape invariance is a powerful symmetry which 
leads to a determination of the spectrum of the system. 
For example, it follows from ^ that the n th energy level 
of the original Hamiltonian iJ_(ao) is given by (the sum 
of the differences in the energy levels) 



En = Y.Riak). 



(5) 
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Shape invariance has been exploited quite successfully 
in various areas of physics (see, e.g., Refs. [TTI - [T7| ). Har- 
monic oscillator is the simplest example of a shape in- 
variant system. In this case, the parameter gq — uj which 
does not change in the sequence of Hamiltonians (namely. 
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FIG. 1: Spectrum for a shape invariant family of Hamiltoni- 
ans. The dashed hnes represent the states that the partner 
Hamihonians do not share. 



ak = /(flfc-i) = flfe-i = w, fc = 1, 2, • • • ), only the poten- 
tial (Hamiltonian) shifts leading to (the shift denoting 
the spacing between energy levels is a constant indepen- 
dent of the level, R{ak) — — lo) 

k 



(6) 
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As one knows, the harmonic oscillator is described by the 
Hermite equation and this analysis shows that shape in- 
variance is the underlying symmetry responsible for its 
solubility. In a similar spirit, we can ask if equations 
describing other orthogonal polynomials may also have 
an underlying shape invariance symmetry responsible for 
their nice properties. In this letter we focus only on the 
Legendre and the associated Legendre equations and an- 
alyze the shape invariance associated with such systems. 
As we will see, the operators associated with the Leg- 
endre equation are not of conventional form and, there- 
fore, the discussion of supersymmetry and shape invari- 
ance generalizes slightly. Since these equations as well as 
their regular solutions, namely, the Legendre and the as- 
sociated Legendre polynomials appear in many areas of 
physics whenever we have a system with spherical sym- 
metry (for example, in the multipole expansions etc) , this 
analysis is of general and current interest. 

The Legendre equation can be written as {x = 
cos 6*, -1 < x < 1) 



dx ^ dx 



n{n + l) Pnix)=0, 



where n = 0, 1, 2, 



For simplicity, we define 



H ^ / 2 1) 

dx dx ' 



so that we can write ([7]) also as 

(i?-n(n + l))P„ = 0, 



(7) 



(8) 



(9) 



From the recursion relations satisfied by the Legendre 
polynomials, we can determine the raising and lowering 
operators for the Legendre polynomials to be 



(x^ -l)-^+nx= {x^ - l)"^+i-^(a;2 
dx dx 

_ A(^2 _ 1) + „^ = _(.^2 _ l)t_i(2,2 

dx dx 



so that we have (for n > 1 and Pq = 1) 
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(10) 



(11) 



and satisfy the equation 

(a„+24 - "')^'n-i = = (1 - a;2) (if - [n - l)n) P„^i. 

(12) 

This equation can also be written equivalently as 
(a|j_]^a„+i — (n — l)^)P„_i = and would generate the 
polynomials P_„_i which are not independent so that 
we will continue with (I12p . In other words, these poly- 
nomials indeed satisfy the Legendre equation. We note 
here that, unlike the harmonic oscillator system, here the 
forms of the raising and lowering operators depend on 
the level leading to a nonuniform (n dependent) spacing 
between the levels as we will see. 

The Legendre polynomials can be iteratively con- 
structed from pT|) and with the help of (|10p we obtain a 
closed form expression of the form (n > 1) 



P„ 



1 



(x^-iy 
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dx 



(13) 

where we have used Pq = 1. This form of the polyno- 
mial directly leads to the recursion relations, but seems 
quite different from the Rodrigues' formula [l^ which 
also gives a closed form expression for these polynomials. 
For low values of n, one can directly check that the two 
expressions are equivalent. To show equivalence in gen- 
eral, let us assume that the two expressions are the same 
up to P„_i so that we can write for n ^ n — 1, also 
as 



Pn-i — 



1 



^(a;2 - 1)'^ 



2"-i((n - 1)!) 
Then, using the identity 
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(14) 
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it can be shown that ([TT|) leads to 
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-i^' ~ 1)- 
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1 d"(a:2 - 1)'^ 



dx 



dx"" ^ 

(15) 

.d"-i(x2 - l)"-i 



dx^' 



2"(n!) 



(16) 
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so that the expression (1131) is indeed equivalent to the Ro- 
drigues' formula, but may be more useful in some cases. 

The supersymmetry of the system is easily seen by 
noting from that if 

a„+2at|V) = i?IV'), (17) 
then for E^Q (or n^Q), 

ata„+2(4l^)) =S(4l^)), (18) 

so that this pair of Hamiltonians can be thought of as 
supersymmetric partner Hamiltonians. However, the dis- 
advantage of this description is that these Hamiltonians 
are not Hermitian. A Hermitian description is not neces- 
sary, but is much more desirable from the point of view 
of supersymmetry and shape invariance and this can be 
done in the following way. 

Let us note from p2)) that we can define a sequence of 
Hermitian Hamiltonians of the form 

H^''-^'^ = (l-a;2)-i (ak+2ai - k^) = H-{k-l)k, (19) 



where fc = 1,2,--- so that ^ H ^ ^{x^ - 1)£. 

Since H does not contain any potential term, the Hamil- 
tonians in ()19|) are shape invariant with only a shift, as 
in the case of the oscillator, but in this case the shift 
depends on the level k. We note here that although the 
left hand side of ([T9ll does not appear to be manifestly 
Hermitian, it really is, following from the identities 



(l-a;2) "a,, = a„_2m(l-a;2)" 



'l-X 



2\-m 



(20) 



where |m| = 0,i,l,-- - These identities are indeed 
quite useful in an algebraic study of this system. In fact, 
using ([20)1 . we can also write i/C'^i) in the manifestly 
Hermitian form 

i/C'-^' = (l-x^)-^ (au+,al^,-k^) {l-x^)-'^. (21) 

The sequence of shape invariant Hamiltonians is now eas- 
ily seen to satisfy 

(1 - x^)-'^ (4+iafc+i - fc") (1 - a^')"^ 

= (1 - ^')"^ (afe+24+2 + If) (1 - x^)--^ 
= ijC^) = i/C^-i) - 2k. (22) 

As we have mentioned before, here the order k in the se- 
quence (level) becomes the relevant parameter and com- 
paring with the harmonic oscillator, we note that the 
shift between two levels now depends on the level. This 
has its origin in the fc-dependent commutation relation 



[afc+i,4+i] = 2fc(l - x^). 



(23) 



which can be contrasted with the universal (level inde- 
pendent) commutator 



[a, a' 



(24) 



for the harmonic oscillator. The n th level oi H = iJ*^"^ 
can now be determined to be 



E„ =2^k = n{n + l), 



(25) 



fc=i 



as expected. This analysis brings out the intimate rela- 
tionship between shape invariance and the solubility of 
the Legendre equation. 

Let us next consider the associated Legendre equation 
which can be written as 



H 



1-x^ 



Pn,m{x) = n{n + 1) 

), (26) 



where m = 0,±1,±2,--- , ±n and H is defined in ([5]). 
This equation is symmetric under m -^r^ —m and so we 
concentrate only on the values m > 0. We note that for 
m — 0, this equation reduces to the Legendre equation in 
0. To analyze the supersymmetry and shape invariance 
associated with the associated Legendre equation, let us 
define the raising and lowering operators 



Al^{l-x') 

„2 



dx 



a-x' 



.... d mx 
(1 -a;^)3— + 

dx (1_2;2)2 



dx 



d , 1 mx 

— {I - x^)^ + ■ 

dx^ {1-x^f 



so that we have the recursion relation 
With these we can now write (see (|26|) ) 



(27) 



(28) 



= H + 



1 - x^ 



= A^mAl + m{m + l), (29) 



where ^ H is defined in 



We note that for any 
{AraAl+mim + l))\,p)=E\i;), (30) 



then, 



(4,^777 + m(m + 1)) (Aim) = E (4,1 V)) , (31) 

so that these can be thought of as supersymmetric part- 
ner Hamiltonians. In fact, it is straightforward to show 
that (see ([29]) ) 



Al^A„, + m{m + I) = A^+iAI^^-^ + (m + l)(m + 2) 

(m -I- 1)2 



1 -a;2 



(32) 
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This shows that at any level of H , namely, En = n{n + 
1), n ^ 0, we can construct a sequence of shape invariant 
Hamiltonians given by 



m = 0,1,2,- •• ,n, (33) 



where the parameter on which the potential depends can 
be identified with = m, a^+i = f{o-m) — m + 1. Fur- 
thermore, this is an example of shape invariance without 
any shift in the Hamiltonians, R{am) ~ 0, so that all 
the (2n + 1) states (including the states with m < 0) 
will have degenerate eigenvalues. Degeneracy in the az- 
imuthal quantum numbers is normally understood as a 
consequence of rotational symmetry in a system. How- 
ever, here we have a new interpretation of this degeneracy 
from the point of view of shape invariance of the system. 

An immediate consequence of this shape invariance 
analysis is that we can relate the associated Legendre 
polynomials, Pn,m, directly to the Legendre polynomials 
through the lowering operators. Constructing the wave 
functions iteratively using (j28p we obtain 



2«(n!) 



(1 



] 2 )^ — 



(34) 



where we have used the definitions in (|27l) as well as the 
Rodrigues' formula ([T6l) . These relations are well known 
in the literature, but shape invariance leads to a simple 
and direct derivation for them. We note here that some- 
times the associated Legendre polynomials (for to > 0) 
are defined with a Condon-Shortley phase (—1)"* which 
can be easily incorporated into our formalism by chang- 
ing the signs of the lowering and the raising operators in 



([27)) . For completeness we note that the associated Leg- 
endre polynomials for to < are proportional to those 
with m > and the exact relation is given by 



(?i — to)! 

Pn,—m{x) — (~1) 7 ; \l P" 

[n + m)\ 



.(x), (35) 



which also follows from (p4|) . 

To conclude we recall that the Legendre and the associ- 
ated Legendre equations are fundamental in the study of 
physical systems with spherical symmetry and their so- 
lutions, the Legendre and the associated Legendre poly- 
nomials, have deep relations in physics. These equations 
have been studied exhaustively in the past. In this letter 
we have described a new way of looking at these equations 
and their solutions from the point of view of an under- 
lying supersymmetry and shape invariance. This can be 
thought of as an operator analysis of the system, much 
like in the case of the harmonic oscillator, which is ob- 
tained by exploiting the underlying supersymmetry and 
shape invariance of the system. This analysis leads di- 
rectly to a closed form expression for the Legendre poly- 
nomials which is equivalent to but is manifestly differ- 
ent from the Rodrigues' formula and which may be use- 
ful in some studies. Degeneracy of the azimuthal quan- 
tum numbers in the case of rotationally invariant sys- 
tems finds a new interpretation from the point of view 
of shape invariance. Other properties such as the rela- 
tion between the associated Legendre polynomials and 
the Legendre polynomials follow directly from such an 
operator analysis. There are, of course, other orthogonal 
polynomials which are also of importance in physics. It 
would be interesting to examine which of these are de- 
scribed by equations with an underlying shape invariance 
symmetry. 
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